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i heating  of  the  ground  by  the  microwave  beam  providing  a source  of  infrared 
power  for  dissipating  the  fog. 

The  study  showed  that,  for  either  caseA very  large  power  densities,  well 
above  the  personal  safety  limit  of  100  W/m^  used  in  the  U.  S. , would  be 
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A Feasibility  Study  of  the  Use  of 
Radiant  Energy  for  Fog  Dispersal 


1.  INTRODUCTION 

. 

For  over  15  years  the  Air  Force  has  been  studying  and  evaluating  various 

methods  of  dispersing  warm  fog.  Many  possibilities  have  been  investigated, 

including  sensible  heat,  helicopter  downwash,  hygroscopic  seeding,  carbon  black, 

and  electric  charging  of  fog  particles.  One  possible  method  that  has  surfaced  in 

recent  years  is  the  use  of  microwave  energy  to  heat  the  water  vapor  and  water 

droplets  in  the  fog  and  thereby  dissipate  it. 

A preliminary  estimate  indicated  that  the  amount  of  microwave  energy  required 

to  vaporize  the  water  droplets  in  a typical  fog  over  an  airport  runway  would  be  at 

least  of  the  order  of  10**  J.  In  addition,  it  appeared  highly  likely  that  the  power 

density  required  to  clear  the  fog  in  a period  of  several  minutes  would  exceed  the 

2 

maximum  power  density  of  100  W/m  that  is  generally  employed  in  this  country  for 
personal  safety.  * At  this  power  level  the  time  of  exposure  recommended  is  6 min 
or  less. 

Despite  these  preliminary  findings,  it  was  felt  that  a more  detailed  study 
would  be  useful  to  help  assess  more  precisely  the  feasibility  and  utility  of  such  a 

(Received  for  publication  5 October  1978) 

1.  Baranski,  S. , and  Czerski,  V.  (1976)  Biological  Effects  of  Microwaves. 

Dowden,  Hutchinson  & Ross,  Inc.,  Stroudsburg,  Pa.,  Chap.  6,  pp  170-173. 
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scheme  ir.  the  light  of  present-day  technology.  Such  a study  would  also  provide  a 
useful  reference  in  the  event  of  an  advance  in  present-day  technology  that  signifi- 
cantly reduces  the  cost  of  microwave  power. 

The  source  of  microwave  energy  was  considered  to  be  ground-based  or  air- 
borne. For  simplicity  a one -dimensional  model,  that  is,  plane-wave  propagation, 
was  assumed.  For  the  ground-based  source,  the  beam  was  taken  as  parallel  to 
the  ground  and  along  the  runway;  it  is  partially  absorbed  by  the  fog,  but  it  does  not 
interact  with  the  ground.  For  the  airborne  microwave  source,  the  beam  was  taken 
as  perpendicular  to  the  ground;  here  the  beam  is  absorbed  both  by  the  fog  and  the 
ground.  The  reflected  microwave  beam  provides  further  heating  of  the  fog,  while 
the  microwave  heating  of  the  ground  provides  an  additional  source  of  heating  of  the 
fog.  The  airborne  method  obviously  provides  a much  better  utilization  of  the 
microwave  energy  than  does  the  ground-based  approach. 


2.  GROUND  BASED  SOURCE:  BEAM  PARALLEL  TO  GROUND  AND 
ALONG  RUNWAY 


For  small  wavelengths  comparable  in  size  to  the  liquid  droplets  in  a fog, 
that  is,  infrared  radiation,  the  scattering  and  absorption  of  the  incoming  radiation 
depends  strongly  upon  the  size  distribution  of  the  fog  droplets.  For  microwave 
frequencies,  however,  the  droplet  size  is  small  compared  to  the  wavelength  and, 
therefore,  the  attenuation  becomes  independent  of  the  size  distribution  and  may  be 

specified  with  respect  to  the  liquid  water  content  of  the  fog. 

o 3 2 

For  a cloud  at  10°C  with  liquid  water  content,  M = 1 gm /m  , Battan  quotes 

the  attenuation  coefficient  K,  dB/km,  for  several  values  of  the  wavelength  A.; 

A = 0.  9 cm,  K = 0.  681 
A = 1.  24  cm,  K = 0.  406 

A = 1.  8 cm,  K = 0.  179. 

o 3 3 

For  inland  fog  at  18°C,  M = 1 gm/m  , Chen  gives: 


A = 0.3  cm,  K = 4.  6 
A = 1 cm,  K = 0.  43 
A = 1.25  cm,  K = 0.25  , 


2.  Battan,  Louis  J.  (1973)  Radar  Observations  of  the  Atmosphere.  University  of 

Chicago  Press,  Chicago  and  London,  Chap.  5,  p.  70. 

3.  Chen,  C.C.  (1975)  Attenuation  of  Electromagnetic  Radiation  by  Haze.  Fo 


For  small  droplets  of  water  at  18°C,  M = 1 gm/m3,  Goldstein4  gives: 

X = 0.2  cm,  K = 7.  14 
X = 0.  5 cm,  K = 1.  65 
X = 1.  25  cm,  K = 0.  280. 

The  attenuation  coefficient  increases  with  a decrease  in  temperature.  At 
X = 1.  25  cm,  the  increase  is  93  percent  in  going  from  18°C  to  0°C.  The  variation 
with  wavelength  is  far  more  sensitive,  a decrease  in  wavelength  by  a factor  of  3 
increasing  the  attenuation  coefficient  by  an  order  of  magnitude.  For  the  purpose 
of  absorption,  a short  wavelength  is  desirable.  However,  in  order  to  generate  a 
fairly  large  amount  of  power,  the  wavelength  cannot  be  much  less  than  1 cm. 
Based  on  the  foregoing  data,  we  shall  use  a value  K = 0.  5 dB/km  per  unit  value  of 
M,  for  wavelengths  near  1 cm,  as  a representative  value  for  fog. 

We  now  consider  the  attenuation  of  a microwave  beam,  parallel  to  the  ground 
and  along  the  runway,  passing  through  a fog  100  m high,  100  m wide,  and  4000  m 

3 

long.  We  shall  assume  a liquid  water  content  M = 0.  5 gm/m  for  which 
K = 0.  25  dB/km.  Since  the  wavelength  is  large  compared  to  the  size  of  the  drop- 
lets, the  absorption  is  large  compared  to  the  scattering,  and  we  take  the  attenua- 
tion of  the  beam  as  due  to  absorption  only.  The  intensity  of  the  transmitted  beam 
relative  to  the  incident  beam  is  given  by 

± 10-kl/10  (1) 

o 

where 

Iq  is  the  incident  energy  flux 
I is  the  transmitted  energy  and 
L is  the  path  length. 

For  a path  L of  4000  m,  corresponding  to  the  length  of  the  runway  and  approach 
zone,  the  attenuation  or  fractional  energy  absorbed  is 


0.2  , 


(2) 


so  that  only  20  percent  of  the  incident  beam  is  absorbed  by  the  fog. 


4.  Goldstein,  H.  (1951)  in  Propagation  of  Short  Radio  Waves  (Kerr,  Editor), 
McGraw-Hill  Book  Co.,  Inc™,  New  York,  Chap.  8,  p.  677. 


With  a liquid  water  content  M = 0.  5 gm/m  , the  mass  of  liquid  water  encoun- 
tered in  a unit  area  of  the  beam  is 


ML  = 2000  gm/m 

For  water  at  20°C,  the  heat  of  vaporization  V is  2450  J/gm;  therefore,  the  energy 
density  G required  to  vaporize  the  fog  is 


which  yields 


G = 2.45  X 107  J/m2  . 


The  temperature  rise  required  to  accommodate  the  additional  water  vapor  is 

about  1°C,  with  a corresponding  heat  input  roughly  equivalent  to  that  required  to 

7 2 

vaporize  the  fog.  Thus,  the  total  energy  density  required  is  about  4.  9 X 10  J/m  . 

If  the  time  allowed  to  clear  the  fog  is  taken  as  10  min,  the  required  power 
density  S is 


4.  9 X 10  n . , 2 

600  82'  000  W/m 


This  is  a very  large  power  density,  well  above  the  safety  limit  in  the  U.  S.  of 
2 

100  W/m  . Assuming  the  fog  is  cleared  to  a depth  of  60  m,  the  cross  sectional 


3 2 

area  normal  to  the  beam  is  6 X 10  m . 
during  this  period  of  10  min  is 


Accordingly,  the  total  power  P required 


P = 6 X 10°  X 82,  000 
= 492  MW  . 

This  represents  a total  energy  requirement  of  about  3 X 1011  J. 

Because  of  air  movement,  a dispersal  time  of  10  min  would  have  to  be  taken 
as  the  maximum  practically  allowable.  From  a realistic  point  of  view,  a system 
capable  of  dispersing  fog  under  most  observed  wind  conditions  would  have  to  be 
able  to  disperse  the  fog  in  a period  of  1 to  2 min.  Thus,  the  time  requirement  for 
the  fog  dispersal  system  considered  here  makes  the  amount  of  power  required 
unrealistically  large  and  incompatible  with  the  maximum  power  density  dictated 
by  personal  safety  regulations. 
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3.  AIRBORNE  SOURCE:  BEAM  PERPENDICULAR  TO  GROUND 


3.1  Analysis  of  Amount  of  Beam  Reflected  by  Ground 

In  this  situation  the  microwave  beam  passes  through  a 100 -m -deep  fog  before 
reaching  the  ground.  Using  Eq.  (1)  we  obtain  a small  absorption  of  about  1 percent. 
The  beam  reaching  the  ground  is  partially  absorbed  by  the  ground,  the  remainder 
being  reflected  back  through  the  fog.  The  surface  of  the  ground  will  heat  up,  the 
temperature  rise  becoming  negligible  when  the  heat  transferred  by  radiation  and 
convection  balances  the  energy  flux  in  the  incident  beam. 

A study  of  the  reflection  characteristics  of  microwave  radiation  shows  a wide 
variation  in  the  amount  absorbed,  depending  on  the  polarization  and  frequency  of 
the  beam,  the  type  of  material,  amount  of  moisture  present,  and  whether  the  sur- 
face is  rough  or  smooth.  However,  in  general,  a large  fraction  of  energy  ranging 
from  50  to  100  percent  is  absorbed.  For  the  purposes  of  this  investigation,  we 
shall  assume  that  100  percent  of  the  energy  is  absorbed.  It  will  be  seen  that  the 
precise  amount  chosen  has  no  qualitative  effect  upon  the  conclusions  drawn  from 
this  study. 


3.2  Determination  of  Temperature  of  Ground 

Although  the  solution  to  the  problem  of  determining  the  temperature  history  of 
a surface  both  radiating  and  subject  to  an  incoming  flux  is  not  given  directly,  the 
problem  may  be  reduced  to  one  of  surface  radiation  only,  for  which  the  solution  is 
known,  by  a suitable  change  of  variable.  We  shall,  as  is  customary  in  problems 
of  this  type,  where  the  temperature  rise  is  not  excessive,  assume  the  linear 
Newtonian  radiation  law  rather  than  the  more  complicated  Stefan-Boltzmann  4th 

5 

power  law  (Carslaw  and  Jaeger  ).  It  may  be  noted  here  that  use  of  the  Newtonian 
law  is  equivalent  to  using  an  average  value  of  the  coefficient  of  convective  heat 
transfer,  H,  which  may  be  too  low  when  the  temperature  rise  is  not  negligible.  It 
will,  therefore,  be  useful  to  obtain  the  effect  of  changes  in  H upon  the  solution. 

We  then  wish  to  solve  the  heat  conduction  equation 


9T 

at 


92T 

° 2 

ax'* 


(5) 


5.  Carslaw,  H.  S. , and  Jaeger,  J.C.  (1959)  Conduction  of  Heat  in  Solids.  Oxford 
University  Press,  London,  Chap.  I,  pp  18-19. 
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where 


PC. 


is  the  coefficient  of  thermal  diffusivity 


p is  density 
C is  specific  heat 

k is  the  coefficient  of  thermal  conductivity 

x is  the  distance  measured  normal  to  the  surface,  and  T is  temperature. 
At  the  surface,  x - 0,  we  have  the  boundary  condition 


9T 

k^T  = H(TS  -Ta)  -F,  x = 0 


(6) 


where  Ta  is 


the  ambient  temperature,  Tg  the  surface  temperature,  and  F U il.c 
incoming  power  flux  density.  We  assume  here  initial  equilibrium,  that  is,  the 
initial  surface  temperature  is  the  same  as  the  ambient. 

We  may  write  Eqs.  (5)  and  (6)  in  the  form 


9v  92v 
9t  = ° 9x2 


(7) 


and 


9v 

9x 


h v 


x = 0 


(8) 


where 


v = T -Ta 


_F 

H 


(9) 


and  h = H/k.  Thus,  a solution  for  v immediately  yields  the  temperature  T since 

F and  H are  known  quantities. 

To  solve  Eqs.  (7)  and  (8)  we  utilize  the  Laplace  transform, 


(10) 
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which,  when  applied  to  Eqs.  (7)  and  (8),  yields 


2- 
q v 


V 


a 


(11) 


and 


= h v , x = 0 , (12) 

dx  ' ’ 

where 

2 / 
q = p/a 

and  V is  the  initial  value  of  v.  Solving  Eq.  (11)  we  obtain 

v = A e + X , (13) 

r 

where  A is  an  arbitrary  constant.  The  positive  exponential  does  not  appear  in 
Eq.  (13)  because  the  temperature  v is  bounded  for  large  x.  To  determine  A,  we 
utilize  the  radiation  condition,  Eq.  (12),  which  gives: 


A - v 

* ~ ' ¥ h + q 


(14) 


and 


v 


V _h_  -qx  V 
P h + q p 


(15) 


Using  a table  of  Laplace  transforms,  we  obtain 


JL  = 
V 


erf 


x 

2(at)1/2 


hx+h2ot  _ 

+ e erfc 


L 2(at) 


1/2 


+ h(at) 


1/2 


(16) 


where  erf  denotes  the  error  function  and  erfc  the  complementary  error  function. 
The  surface  temperature  v , corresponding  to  x = 0,  is  then  given  by 


s*1  a ^ erfc  (h(at)','/^l 


(17) 
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For  the  case  of  no  convective  heat  transfer  (H  = 0),  the  temperature  history  is 

g 

given  by  (Carslaw  and  Jaeger  ) 


For  the  purposes  of  this  study,  we  have  chosen  earth  and  concrete  as  typical 
surface  materials  for  which  the  thermal  data  are  (average  representative  values 
have  been  utilized): 

earth: 

a - 0.  002  cm2/sec 
k = 0.  002  cal/cm  sec  °C 
concrete: 

2 

a = 0.  01  cm  /sec 
k = 0.  005  cal/cm  sec  °C 

Since,  as  will  be  evident  from  the  calculations,  the  solution  is  more  sensitive  to 

n 

H than  to  oork,  a range  of  values  has  been  chosen  for  H.  Schneider  gives  as  a 

typical  value  for  earth,  H = 5 Btu/ft2  hr  °F  = 6.  8 X 10'4  cal/cm2  sec  °C,  while 

for  concrete  and  sandstone  a representative  value  is  0.  65  Btu/ft2  hr  °F  = 8.8  X 10‘5 
2 o 8 

cal/cm  sec  C (Brown  and  Marco  ).  When  using  the  linear  radiation  law,  Eq.  (6), 
the  coefficient  H,  as  indicated  by  Carslaw,  is  an  overall  value  including  the  effects 

of  radiation  and  convection.  Since  these  effects  are,  in  general,  of  the  same 

-3  2 

order  of  magnitude,  the  coefficient  H may  attain  values  as  high  as  10  cal/cm 

o -4  -3 

sec  C.  We  have,  therefore,  chosen  a range  of  values  of  H from  10  to  10 

cal/cm2  sec  °C. 

The  temperature  histories  for  earth  and  concrete  are  given  in  Figures  1 and  2 
where  the  temperature  rise  AT  = T - Ta  is  plotted  against  time  for  several  values 
of  H.  An  incident  flux,  F,  of  1000  W/m^  was  chosen  for  the  calculations.  Since 
the  temperature  rise  is  proportional  to  F,  the  solution  for  any  other  value  of  F can 
easily  be  obtained  from  these  results.  The  results  for  earth  and  concrete  are  very 
close,  so  very  little  difference  in  the  temperature  rise  may  be  expected  because 


6.  Carslaw,  H.  S. , and  Jaeger,  J.C.  (1959)  Conduction  of  Heat  in  Solids.  Oxford 

University  Press,  London,  Chap.  II,  p.  75. 

7.  Schneider,  P.J.  (1955)  Conduction  Heat  Transfer.  Addison-Wesley  Publishing 

Co.,  Inc.,  Cambridge,  MA,  Chap.  10,  p.  265. 

8.  Brown,  A.  I.,  and  Marco,  S.  M.  (1958)  Introduction  to  Heat  Transfer.  McGraw- 

Hill  Book  Co.,  Inc.,  New  York,  Chap.  11,  p.  194. 
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of  variations  in  the  thermal  properties  of  the  surface.  Changes  in  H have  a much 
more  significant  effect  upon  the  temperature  rise.  For  example,  in  a time  interval 
of  1 hr.  Figure  1 (earth)  indicates  the  temperature  rise  increases  from  16°C  to 
32°C  when  H decreases  from  10  2 to  10 

3.3  Interaction  of  Radiant  Energy  With  Fog 

In  standard  treatments  of  infrared  technology  it  has  been  customary  to  assume 

that  absorption  is  small  compared  to  scattering  during  the  interaction  of  infrared 

9 10 

radiation  with  water  droplets  in  the  atmosphere.  ’ Hence,  only  scattering  has 

been  considered  in  determining  the  attenuation,  by  water  droplets,  of  an  infrared 

beam  passing  through  the  atmosphere.  Recent  studies,  however,  have  indicated 

that  absorption  of  infrared  energy  by  water  droplets  in  the  atmosphere  may  be  a 

11  12 

significant  factor  in  the  attenuation  of  an  infrared  beam.  ’ 

Hodges12  shows  that  from  20  to  50  percent  of  the  total  attenuation  of  the  beam 

is  due  to  absorption,  the  higher  amounts  occurring  at  wavelengths  beyond  10  p. 

13 

Since  most  fogs  are  opaque  to  infrared  radiation,  we  shall  assume,  for  the 

purposes  of  this  study,  that  the  total  attenuation  is  100  percent  and  that  40  percent 

of  the  incident  beam  is  absorbed  by  the  fog. 

As  indicated  previously,  the  maximum  safe  power  density  in  the  U.  S.  is 
o 

100  W/m  . For  this  incoming  flux,  the  amount  absorbed  by  the  earth  is,  assum- 

2 

ing  100  percent  absorption,  also  100  W/m  (for  convenience,  we  neglect  here  the 

small  fraction  of  the  incident  microwave  beam  absorbed  by  the  fog).  Taking 
-4  2 

H = 5 X 10  cal/cm  sec  as  a representative  value  for  the  convective  heat  transfer 
coefficient,  the  temperature  rise  for  earth  in  a period  of  1 hr  is,  from  Figure  1, 

AT  = 23°C.  For  a flux  of  100  W/m2,  the  temperature  rise  is  1/10  of  this  amount, 
or  2.3°C.  At  an  ambient  temperature  of  18°C  = 291°K  the  corresponding  surface 
temperature  is  Tg  = 293. 3°K.  Since  a fog  of  high  liquid  water  content  and  a depth 


9.  Kruse,  P.  W.  , McGlaughlin,  P.D.,  and  McQuistan,  R.  B.  (1962)  Elements  of 
Infrared  Technology,  John  Wiley  and  Sons,  Inc.,  New  York,  Chap.  5, 

p.  186.  

10.  Plass,  Gilbert  N.  ( 1965)  Handbook  of  Military  Infrared  Technology  (Wolfe, 

Editor),  Superintendent  of  Documents,  U.S.  Government  Printing  Office, 
Washington,  D.C.,  Chap.  6,  p.  206. 

11.  Carlon,  H.  R.  (1970)  Infrared  emission  by  fine  water  aerosols  and  fogs. 

Applied  Optics  9(No.  9):2000-2006. 

12.  Hodges,  J.A.  (1972)  Aerosol  extinction  contribution  to  atmospheric  attenuation 

in  infrared  wavelengths.  Applied  Optics  lA(No.  10):2304-2310. 

13.  Chen,  C.  C.  (1975)  Attenuation  of  Electromagnetic  Radiation  by  Haze  Fog, 

Clouds  and  Rain.  The  Rand  Corporation,  Santa  Monica,  CA,  R-1094-RR, 

p.  18. 
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of  100  meters  or  more  is  opaque  to  infrared  radiation,  we  may  consider  it  as  a 
radiating  blackbody.  The  net  power  density  flux  Sr  radiating  from  the  surface  to 

the  fog  is  then 


Sr  = ° K4'Ta4 


where 


5.  67  X 10'8  W/m2  °K4  is  the  Stefan -Boltzmann  constant. 


We  then  obtain 


S = 5.  67  X 10'8  [<293.  3)4  - (291)4] 
r 


= 13.  1 W/m2  . 


The  flux  absorbed  by  the  fog,  ASr>  assuming  40  percent  absorption,  is  then 


ASr  = 0.4(13.1) 


= 5.  24  W/m2  . 

For  a fog  depth  of  100  m and  a liquid  water  content  of  0.5  gm/m  , the  heat  G 
required  to  vaporize  the  liquid  droplets  in  a unit  area  of  the  fog,  as  shown  in 
Section  2,  is 

G = VML  = 1.  22  X 105  J/m2  . 

As  in  the  previous  case,  the  amount  of  heat  required  to  raise  the  temperature  of 
the  air  sufficiently  to  accommodate  the  additional  water  vapor  is  about  the  same  as^ 
G.  Thus,  the  total  required  energy  density  is  twice  this  value,  or  2.45  X 10  J/m  . 
The  time  At  required  to  vaporize  the  fog  by  infrared  radiation  is,  therefore. 


2.45  X 10 


= 0.47  X 10  sec 
= 13  hr  . 


Since  the  absorption  for  microwave  energy  is  only  1 percent,  the  direct  effect  of 

2 

microwave  power  at  a density  of  100  W/m  upon  the  fog  dispersal  may  be  neglected. 

We  then  see  that  the  time  required  to  vaporize  the  fog,  when  limited  by  a power 
2 

density  of  100  W/m  , is  unrealistically  long.  Even  if  the  amount  absorbed  by  the 
fog  were  assumed  to  be  100  percent  rather  than  40  percent,  the  time  required 
would  still  be  over  5 hr. 

2 

At  a power  density  of  100  W/m  , the  total  power  required  for  a runway 
4000  m long  by  100  m wide  is  40  MW,  an  extremely  large  amount  of  power  for 
airborne  equipment.  Thus,  even  where  limited  by  the  maximum  power  density 
imposed  by  personal  safety  factors,  the  use  of  microwave  power  by  airborne 
equipment  to  clear  a fog  is  unfeasible  with  regard  to  the  amount  of  power  required. 


4.  SUMMARY  AND  CONCLUSIONS 

A feasibility  study  has  been  made  of  the  possible  use  of  radiant  energy  to  clear 
a fog  over  an  airport.  The  source  of  microwave  power  was  considered  either 
ground-based  or  airborne.  For  the  ground-based  case,  the  microwave  beam  was 
taken  as  parallel  to  the  ground  and  interacting  only  with  the  fog.  For  the  airborne 
case,  the  beam  was  taken  as  perpendicular  to  the  ground,  the  heating  of  the  ground 
by  the  microwave  beam  providing  a source  of  infrared  power  for  dissipating  the 
fog. 

The  study  showed  that,  for  either  case,  unrealistically  large  amounts  of  power, 
well  above  the  personal  safety  limits  prevailing  in  the  U.  S. , were  required  to  dissi- 
pate the  fog  in  a reasonable  length  of  time,  that  is,  about  10  min.  To  conform  to 

2 

the  maximum  power  density  determined  by  safety  requirements,  100  W/m  , 
imposes  unrealistic  times,  that  is,  many  hours,  to  dissipate  the  fog.  In  addition, 
because  of  the  high  cost  of  electrical  energy,  the  large  amount  of  energy  required 
to  dissipate  a typical  airport  fog  (about  3 X 10**  J)  makes  the  scheme  prohibitively 
expensive. 
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